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We introduce the eikonal approximation to study the effect of the large-scale motion of cosmic 
fluids on their small-scale evolution. This approach consists in collecting the impact of the long- 
wavelength displacement field into a single or finite number of random variables, whose statistical 
properties can be computed from the initial conditions. For a single dark matter fluid, we show 
that we can recover the nonlinear propagators of renormalized perturbation theory. These are 
obtained with no need to assume that the displacement field follows the linear theory. Then we 
extend the eikonal approximation to many fluids. In particular, we study the case of two non- 
relativistic components and we derive their resummed propagators in the presence of isodensity 
modes. Unlike the adiabatic case, where only the phase of small-scale modes is affected by the 
large-scale advection field, the isodensity modes change also the amplitude on small scales. We 
explicitly solve the case of cold dark matter-baryon mixing and find that the isodensity modes 
induce only very small corrections to the resummed propagators. 

PACS numbers: 



I. INTRODUCTION 

The development of wide-field surveys has been trig- 
gering renew interest in the implementation of pertur- 
bation techniques for the computation of the statistical 
properties of large-scale structures. Several approaches 
have been proposed to significantly extend the standard 
perturbation theory (PT) methods (see pQ). A particu- 
larly interesting approach is the so-called renormalized 
perturbation theory (RPT), pioneered by Crocce and 
Scoccimarro [2H3]. This method relies on the use of 
the 2-point propagator as a measure of the memory of 
the initial conditions. This appears as the fundamental 
building block from which perturbation theory can be 
re-constructed and allows to take into account nonlin- 
earities from very small scales, reducing their impact in 
the neglected terms of the perturbative expansion. This 
idea was later extended in [5J [5] with the introduction of 
multi-point propagators. 

The key result of RPT is that in the high-fc limit the 
propagators can be computed exactly, by summing up an 
infinite subset of contributions in the standard perturba- 
tion theory expansion. However, this result has been 
proved only for a single pressureless fluid - describing 
cold dark matter (CDM) - using a technique that seems 
difficult to extend to more complex scenarios. Thus, 
there is no systematic way to implement the RPT ap- 
proach when the content of the cosmic fluid is richer - 
e.g. when it includes various matter components, non- 
relativistic neutrinos, or even modification of gravity - 
and its application range has been so far limited to simple 
cosmological models (see however [7]). Note that other 
approaches, such as the so-called time renormalization 
group approach [8., do not suffer from such limitations. 

On the other hand, the resummation of the propaga- 
tors can be obtained by a more direct technique than that 



originally introduced in [2H1]. As mentioned in [5] and 
explicitly used in [6l [10] , in the high-fc limit it is possible 
to resum the same class of contributions making use of 
a single or a finite number of random variables, which 
describe the effect of the long- wavelength fluctuations on 
smaller scales. This has been called the a-method [6]. 
Here we will explicitly present how to compute nonlinear 
propagators in this framework. Moreover, we will show 
that this method can be employed to extend the RPT 
approach to arbitrarily complicated cosmologies. Bor- 
rowing the terminology from quantum field theory, where 
similar techniques are used (see e.g. [HJ[H]), we propose 
to dub this method the eikonal approximation. 

The fluid content of the universe is richer than a simple 
single-dark matter component. In practice we know that 
the properties of the large-scale structure of the universe 
can be significantly affected by the presence of a sub- 
dominant species. This is the case of baryons, which at 
high enough redshift behave very differently from CDM. 
Indeed, the net result of this different behavior is the 
existence of the baryonic oscillations. 

From a theoretical point of view, the CDM-baryon sys- 
tem is very appealing: After decoupling both components 
are pressureless fluids (at least above the baryonic Jeans 
scale) and thus follow geodesic motion [115]. However, the 
matter fluid as a whole cannot be described as a pres- 
sureless effective fluid. The reason is that the CDM and 
baryon fluids are moving at different velocities - see [13] 
for a study of some of the consequences of this different 
behavior and their possible observational implications. 
Such a velocity dispersion induces an effective anisotropic 
pressure in the total fluid, modifying its equations of mo- 
tion. Thus, one is forced to study the system of coupled 
equations for the CDM and baryons, as previously done 
in [7]. In particular, this is the system that we will ex- 
plore with the help of the eikonal approximation. 
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The plan of the paper is the following. In Section [TT] 
we review the basic concepts of the RPT approach for a 
single CDM fluid and we discuss the eikonal approxima- 
tion in this context. In Section IlIII we extend this dis- 
cussion to the multi-fluid case. In particular, we derive 
the evolution equations for several gravitationally cou- 
pled pressureless fluids, we describe the various modes 
that appear in this case, and we present how they can 
be incorporated in the eikonal approximation. Finally, 
in Section [lV| we illustrate our concepts in the case of the 
standard CDM-baryon mixing. 

II. SINGLE FLUID 

In this section we review the basic concepts of the RPT 
approach with a single perfect fluid, developed in [HH]. 
Moreover, we rederive the procedure to resum the non- 
linear propagators for cosmic fluids using the eikonal ap- 
proximation. 

A. Equations of motion 

We assume the universe to be filled by one pressure- 
less fluid. We denote its density by p, and the density 
contrast by 5 = p/p — 1, where J> is the average energy 
density. The continuity equation then reads 

p+l ((! + *)«% =0, (1) 



where u l is the i-component of the peculiar velocity field 
of the fluid and a comma denotes the partial derivative. 
The Euler equation is 

®-y + Hu l + -uiu l 3 =--<j>,i , (2) 
at a u a 

where H is the Hubble rate, H = dlna/dt, and cj) is 
the gravitational potential. Since we are only interested 
in the dynamics on sub-horizon scales, <j) is the usual 
Newtonian potential, satisfying the Poisson equation 

Acj) = A-kG a 2 p5 . (3) 

We also ignore small-scale shell crossings in the fluid. 
Then, since the gravitational force is potential the fluid 
velocity remains potential at all orders in the perturba- 
tions. Thus, it can be entirely described by the dimen- 
sionless velocity divergence, defined by 

s = ^- < 4 > 

By using the following convention for the Fourier 
modes, 

/(k) 55 / l^ 6 "^ > (B) 

the equations of motion can then be rewritten in Fourier 
space as 



1__9 

h at 



<y(k) + fl(k) = -a(k 1 ,k 2 )d(ki)5(k 2 ) , 



(6) 
(7) 



where f2 m is the reduced matter density and 



B. RPT formulation 



a(ki,k 2 ) = 
/3(ki,k 2 ) = 



(ki + k 2 )-ki 

kl ' 
(k 1 + k 2 ) 2 k r k 2 



(8) 
(9) 



On the right-hand side of eqs. ([6| and integra- 
tion over repeated wave modes and a Dirac function 
5n(k — ki — k 2 ) is implied. Note that these equations 
are valid irrespectively of the dark energy equation of 
state or curvature term. 



In order to recast these equations in RPT form let 
us first discuss their linear solutions. At linear order, 
the coupling terms in the right-hand side of eqs. ^ and 
([7]) are absent. We are then left with the usual linear 
solutions of a pressureless fluid, i.e. 



<J(x, t) = D+ (t)6+ (x) + £>_ (t)6- (x) 



(10) 



where D+(t) and D—(t) correspond to the growing and 
decaying modes, respectively. The corresponding expres- 
sion for the dimensionless velocity divergence is 

0(x,t) = -f + {t)D + (t)5+(x) - /_(t)Z?_(t)<y_(x) , (11) 

where / + and /_ are the growth rates, defined as f± = 
dln.D ± /dliia. 
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As introduced in [TJ] , it is convenient to define the 
duplet 



(12) 



where is the reduced velocity contrast defined as 

e(x,t) = -0(x,t)// + (t) 

f (t) (13) 
= D + (t)5 + (x) + g^D_(t)*_(x), 

in such a way that the linear growing mode of O is the 
same as that of S. It is then convenient to rewrite the evo- 
lution equation using rj as time variable, defined through 



D+drj = dD^ 



(14) 



With this definition the equations of motion ^ and ^ 
can be recast as 



d 



dr) 
where 



*„(k) + n ab V b (k) = 7 ahc (k, ki, k 2 )* 6 (ki)* c (k 2 





3 n m 3 n,„ 
'27+ *T+ 



-1 



1 



(15) 



(16) 



and the non-zero elements of the coupling matrix ^f a bc are 

a(k 2 ,ki) 



7n 2 (k,ki,k 2 
7i 2 i(k,ki,k 2 



a(ki,k 2 ) 



(17) 



7 222 (k,ki,k 2 ) =/3(ki,k 2 ) 



From cqs. (10) and (13), the growing and decaying 



solutions are proportional to 

4+) = (l,l) T and «<-> = (!, /_//+)' 



(18) 



respectively. It has been widely stressed that /-//+ is 
very weakly dependent on the background. Indeed, it 
departs little from the value it takes in an Einstein-de 
Sitter universe (EdS), i.e. f2 m = 1, where /-//+ = —3/2 

The solutions of the linear equations of motion, ob- 



tained from dropping the right-hand side of (15), can 



be formally written in terms of the linear propagator 
gabiViV')- This is such that 



and 



gab{V,V) = S ab , (19) 



9ab{v, Vo) + ^ac(v)9cb(ll, %) = . (20) 



It can be built from a complete set of independent solu- 
tions of the evolution equation. For a single fluid we can 



use ui + ^ and Ua ' defined in eq. (1 1 8h . We then have 



,(-) 



tofMoJ^^W^fe), (21) 



vf^fkl 



*i 2) (k) 



g a b{g,go) 



-O *b(k) 



ce(v',V0)^O *e(kl) 




9df(v', VoF^O *f(k 



FIG. 1: Diagrammatic representation of the series expansion 
of ^ a (k) up to fourth order in the initial conditions denoted 
here by $ a (k). Time increases along each segment according 
to the arrow and each segment bears a factor g c d(gf — rji) 
if r]i is the initial time and gj is the final time. At each 
initial point and each vertex point there is a sum over the 
component indices; a sum over the incoming wave modes is 
also implicit and, finally, the time coordinate of the vertex 
points is integrated from 770 to the final time g according to 
the time ordering of each diagram. 



where the coefficients are chosen such that 



£ «a Q) fa) <r = 



(22) 



For an EdS background the explicit form of g a b reads 
e v-vo ( 3 2 \ e-ffo-'w) [2-2 



g a b{g,vo) 



3 2 



-3 3 



The linear propagator is useful to formally write the 
solution for ty a in integral form. Indeed, using eqs. (12) 
and (13), the equations of motion ^ and Q can be 
written as [2] 

*a(k, rj) = g ab (r ll ry )*b(k, 770) 

dv'gabiv, v'hbde(K ki, k 2 )* d (ki, f/)*e(k 2 , 77') . 

(24) 

As illustrated in Fig. [I] this equation has a diagrammatic 
representation in the RPT context [2J. 

Another important quantity introduced in the RPT 
approach is the nonlinear multi-point propagator. More 
precisely, the (n+ l)-point propagator b is defined 

by 



9* i , 1 (ki ) r ?0 )...^ 6n (k„,77o) 

n 

s o(k - k *) r ibl...&„( k i' • • ■ > k «; go) 



(25) 



Propagators represent the way the ^ a 's respond to an 
infinitesimal change of the modes at an earlier time and 
they are important in the construction of multi-point 
spectra 0(3]. 
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In the large-fc limit (to be better specified below) these 
propagators enjoy a remarkable property. Indeed, in [3] it 
has been shown that in this limit and for Gaussian initial 
conditions, the nonlinear 2-point propagator G a b 
has a simple expression 



= r« 

1 ab 



G a b(k; f], r) ) = g a b(v, Vo) exp ( 



-k 2 a 2 A 



J ) 2 /2) 



(26) 

where cr^ is the variance of the initial displacement field. 
Note that the linear propagator <? a f, is simply the tree 
level analog of G a b- This result has been generalized to 
(n + l)-point propagators with n > 2 in [5], where it has 
been shown that 



= r 



(n)— tree 
ab\ ...b n 



exp I 



-fc 2 ^(e" - , (27) 



where r(™) _troc is the corresponding propagator com- 
puted at tree level. 

The exponentiation in eqs. (26) and (27) has been ob- 



tained in [21 [5] by summing up an infinite number of 
diagrams thought to dominate in the large-fc limit. In 
order to identify which diagrams dominate in this limit, 
the concept of principal line and its generalization for the 
(n + l)-point propagators, the principal tree, have been 
introduced. In 3] it has been shown that each diagram 
contributing to the nonlinear propagator G a b(k; r\, 770) al- 
ways contains a unique line that goes from some time rj 
(symbolized by the vertical dotted line) to a final time 
77. To this line may be attached loops containing power 
spectra evaluated at an initial time rji n . This is illus- 
trated in Fig. [2j upper panel. The principal line is the 
unique way to go from 770 to 77 without crossing an initial 
power spectrum ®, thus moving always in the direction 
of increasing time. Similarly, for each diagram contribut- 



ing to r 



(«) 

ab\ ...b r . 



there always exists a unique tree with n 



branches, the principal tree, that joins 770 to 77 (see bot- 
tom diagram of Fig. [2]) [5] . 

We can now specify under which assumption the rela- 
tions (26) and (27) have been derived. These are: 



• The multi-point propagators are dominated by 
those diagrams in which every loop is directly con- 
nected to the principal tree. 

• The diagrams are computed and summed up in the 
limit where the incoming wave modes qi are soft, 
i.e. qi -C fc. 

As we will show below, the eikonal approximation corre- 
sponds exactly to the last assumption. It can incorporate 
the first one if necessary. 



Resumming the 2-point propagator with the 
eikonal approximation 



In [6] it has been shown that eqs. (26) and (27) can 
be obtained irrespectively of the diagrammatic repre- 
sentations and of the nature of the initial conditions. 



Gab(k) - 



ri 3 b L(ki,k 2 ,k 3 ) k 




FIG. 2: 



andrg^k! 



Example of diagrams contributing to G a b(k) (top) 
k2,ka) (bottom). The dominant contribution 
after resuming all possible configurations is expected to come 
from those diagrams where all loops are directly connected 
to the principal line (top) or principal tree (bottom). The 
principal line and tree are drawn with a thick solid line. A 
symbol <g> denotes a power spectrum evaluated at initial time 
rji n . The dominant loops are those drawn by dashed lines, 
while the sub-dominant loops are those in dotted lines. 



Indeed, the nonlinear fluid equations contain nonlinear 
terms that couple short and long-wavelength modes. The 
eikonal approximation corresponds to study the effect of 
very long-wavelength modes q on the dynamics of a given 
short-wavelength mode fc, in the limit of q <C fc. In this 
limit, space variations of the long- wavelength modes are 
tiny with respect to the mode fc, and the long modes 
can be treated as an external random background. If wc 
neglect the mode couplings between short scales, the non- 
linear fluid equations can be rewritten as linear equations 
embedded in an external random medium. 

Let us be more explicit here. Coupling terms are given 
by a convolution of fields taken at wave modes ki and 
k2 such that k = ki + k2. These nonlinear terms can 
be split into two different contributions: the one coming 
from coupling two modes of very different amplitudes, 
fci "C fc2 or &2 <C fci, and the one coming from coupling 
two modes of comparable amplitudes. In the first case, 
the small wave modes ought to be much smaller than k 
itself. Let us denote these small modes by q. In the limit 
of q <C fc, the equations of motion ( 15 ) can be rewritten 
as 



JU a (k) + fU* h (k) = H a6 (k)* b (k) 
+ [ 7a6c (k,ki,k2)* 6 (k 1 )* c (k 2 )] w , 



with 



S ab (k,7i) = 2 / d 3 q7 Q6c (k,k,q)* c (q,77) 
Js 



(29) 
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The key point is that in eq. ( 29 1 the domain of integra- 



tion is restricted to the soft momenta, for which q <C k. 



Conversely, on the right-hand side of eq. ( 28 ) the convo- 



lution is done excluding the soft domain, i.e. it is over 
hard modes or modes of comparable size. 

In the limit of separation of scales, 2 a b is a random 
quantity which depends on the initial conditions. Us- 
ing eqs. pi and ffify, for q <C k the leading expres- 
sion of the coupling matrix is obtained with the fol- 
lowing limit values a(q, k) (q • k)/g 2 , a(k, q) w 0, 
/3(q,k) = /3(k,q) « (q- k)/(2g 2 ). Thus, labc in eq. ph 



simplifies and S t becomes proportional to the identity, 
with 



E ab (k,i]) = S(k,?7) S ab 
S(k,T7) = 



,3 k-q (30) 
d J q —j- 6(q, rj) . 
s 1 



Note that only the velocity field 9 (and not the density 
field 5) contributes to S a h. Furthermore, as 0(x, 77) is 
real 0(— q) = 6*(q) and thus S is purely imaginary. 



In eq. ( 28 ) we have reabsorbed the effect of the non- 



linear coupling with long-wavelength modes in the linear 
term Safj^t,. The solution to this equation can be given 
in terms of the resummed propagator £ a &(k, 77, rj') [21] 
satisfying the equation 

d_ 

and reads 



S(k, rj) ) U(k, V, rf) + n ac (r/)U(K rj, rf) = 0, 

(31) 



dv'U(ihv') [76rf e (k,k 1 ,k 2 )* d (ki,?7')* e (k2,r]')]« • 

(32) 



where in the last line the convolution is done on the hard 
domain %. 



In the case of a single fluid, as discussed here, eq. (31 ) 



can be easily solved. Taking into account the boundary 
condition £ a &(k, 77, rf) — 5 a b, one obtains 

£ a6 (k, 77,770) = g ab (v,Vo) cx P (^J drf 3(k, rj')j • (33) 

The argument of the exponential is the time integral of 
the velocity projected along the direction k, i.e. the dis- 
placement component along k. Note that in their original 
calculation, Crocce and Scoccimarro assumed that the in- 
coming modes in the soft (i.e. large-scale) lines were in 
the linear and growing regime. Here we need not make 
this assumption. Equation (33) is valid irrespectively of 



the fact that the incoming modes in H are in the growing 
mode or not. 



There is another important aspect of eq. (33). Since 



H(k, 77) is a purely imaginary number, the soft modes 
change only the phase of the small-scale modes but not 
their amplitude. Such an effect will then have no impact 



on the equal-time power spectra. However, it has some 
impact on the amplitude of the propagators. Indeed, the 
phase change inevitably damps the correlation between 
modes at different times. This effect is at the heart of the 
regularization scheme used by approaches such as RPT. 

To illustrate this last point, let us see how one can 
recover eq. ( 26 1 using the solution ( 32 ) and the resummed 



propagator ( 33 ) derived with the eikonal approximation. 
Deriving eq. (32) with respect to an initial field ^(k, 770) 



as in eq. (25 1, and taking the ensemble average one finds 



Gab(k,ri,ri ) = (£ ab (k, 77, rj ))z 



(34) 



The nonlinear 2-point propagator G a b is given by the 
ensemble average of £ &(k, 77, 770) over the realizations of 
3(k). In general, the expression of the nonlinear propa- 
gator introduces the cumulant generating functions of 3. 
Indeed, using (33) eq. (34) yields 



G ab {k,ri,r] ) = g ab (n,rj ) exp V I , (35) 

where c p is the p-order cumulant of the field 
d?7'S(k, 77') and for symmetry reasons the sum is re- 
stricted to even values of p, thus ensuring that the non- 
linear propagators are real. 

For Gaussian initial conditions and assuming that at 
late time the long-wavelength S is in the linear growing 
mode, cumulants with p > 2 in eq. ( 35 1 vanish and we 



are left with only C2. This is given by 

c 2 (k)= /'dr/d77"(S(k,77')S(k,77")) . (36) 
J no 

Then, exploiting the time dependence of the linear grow- 
ing mode, O oc D + = e''~'" n , and using eq. (30), we have 

5(M) = D+( V ) f d 3 q k ^9(q,77 in ) . (37) 
Js 1 



Plugging this expression in eq. ( 36 1 , we can express C2 in 



terms of the initial power spectrum P; n (g), defined by 

(6(q, 77 in )8(q', »fcO) = <5 D (q + q')P in (q) ■ (38) 
Indeed, we have 



c 2 (k) = -k 2 al(e r >- r >"> - e" - 71 '-)' 



(39) 



where a\ gives the variance of the displacement field de- 
fined as [1 13] 



3.// iq T 



(40) 



At this stage a\ depends on the domain of integration 
and hence on k. The standard RPT results are obtained 
by taking the value of a\ in the large-fc limit. We will 
comment on this assumption in the conclusion. Then, 
setting here and in the following ?7i n = for convenience, 
from eq. ( 35 1 we recover eq. (|26|) , 



G ab (k, V, Vo) = 9ab(v, Vo) exp (-fc 2 a 2 (e" - e"«) 2 /2) 



(41) 



G 



D. Higher-order propagators 

Although the focus of this paper is on the 2-point prop- 
agator, let us comment on the use of the eikonal approx- 
imation, in particular of eq. (32), in investigating the 



resummation of higher-order propagators. 

The computation of the nonlinear 3-point propagator 
proceeds by replacing ^ and fy e in the second line of 
eq. (32) by the linear solution given by the first line of 
this equation. Deriving twice with respect to the initial 
field yields 



d 2 * a (k,77) 



•mt, i hi, i , / dr ?'£ad(k; V, v') 
3# 6 (ki,7?o)d# c (k2,i7o) J n 

x [7 de/ (k,ki,k 2 )C e b(ki;ry',?7o)C/c(k2;f7 / ,??o)] w 



(42) 



This is the same formal expression as for the naked the- 
ory except that here the convolution is restricted to the 
hard-mode domain T-L. Note that the coupling vertex be- 
tween modes with hard momenta in the second line is 
not affected by the use of the eikonal approximation: ft 
is identical to the one of the naked theory. Moreover, it is 
remarkable to see that, using the form given by eq. (33), 



the exponential terms factor out of the time integral and 
their arguments sum up to give 



9*i(ki,%)a* c (k a) »jo) 



exp ( [ dry' H(k, rf' 



x / di 1 / g ad (ri,i 1 l )[j def (k,k u k2)geb{v',Vo)9fc{v',Vo)] H - 
J no 

(43) 

Finally, taking the ensemble average and using the def- 
inition of multi-point propagators, eq. (25), one obtains 
in the Gaussian case 



-,(2)-tree. 



\rj im ) exp(~k 2 a 2 d (e^-e^) 2 /2) . 

(44) 

This result can be generalized to propagators of any 
higher order. The formal expressions of the resummed 
trees computed in the eikonal approximation are ob- 
tained from those computed in the naked theory by sim- 
ply changing the propagators from g a b to £ a f,. Then, for 
each pair of merging branches with equal initial time, 



one can factor out the phase similarly to what is done 
when going from eq. ( 42 ) to ([43]) . Finally, this leaves an 



rail factor exp ( dr/' 5(k, ry')^ , which can be factor- 
ized out, recovering eq. (27). The eikonal approximation 



explicitly shows how the results [5J |5] can be recovered 
and generalized to any time-dependent large-scale wave 
mode. 



III. MULTI-FLUIDS 

In this section we explore the case where the universe 
is filled with several non-interacting pressureless fluids 
and show how the eikonal approximation can be imple- 
mented in this case. Due to the gravitational coupling 
and the expansion, at late time such a system becomes 
indistinguishable from a single- fluid component. How- 
ever, during its evolution it can behave very differently 
from a single perfect fluid depending on the initial con- 
ditions. 



A. The equations of motion 

Denoting each fluid by a subscript a, the continuity 
equation reads, for each fluid, 



while the Euler equation reads 

d 1 1 

— < + Hul + -uiu l aj = --4 



(45) 



(46) 



The Poisson equation ([3|, where now S is the density 
contrast of the total fluid energy density, i.e. 



Pm = ^2 Pa = (1 + S m )p n 



(47) 



allows to close the system. This introduces couplings 
between the fluids. 

In Fourier space, the equations of motion become now 



Hdt 



<Uk) + Q (k) = 



Hdt 



'a(k) 



1 dln(a 2 7T) 



H 



dt 



0a(k) + -ft m i m (k) 



-a(ki,k 2 )0 a (ki.)5 a (k2) , 
-/3(k 1 ,k 2 )^(k 1 )^(k 2 ) , 



(48) 
(49) 



r 



where 9 a is the dimensionless divergence of the velocity 
field of the fluid a and J7 m is the reduced total density of 
the pressureless fluids. The coupling between the fluids is 



only due to the term <5 m appearing in the Euler equation. 

Before studying these equations let us discuss the equa- 
tions for the total fluid. As we are describing here a col- 
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lection of pressureless particles, it is tempting to write 
down the equations of motion for the total fluid. The 
continuity equation is simply identical to eq. ([lj. The 
total fluid velocity u l is defined by 

= (50) 

a 

where f a = p a /pm, and its evolution equation reads 

1 , 



d 1 f 

— u l + Hu l + -v?u\ = 4> 4 

at a ' J a ap„ 



where and is the velocity dispersion of the mean fluid, 
given by 



(52) 



Thus, due to the multi-fluid nature of the system, the Eu- 
ler equation contains an anisotropic stress term. One can 
write down an equation of motion for this term, but this 
will involve higher moments of the fluid distribution and 
so on. Thus, the complete description of the total fluid 
at nonlinear order requires an infinite hierarchy of equa- 
tions in the moments of the fluid. Another consequence 
of this expression is that, even though the velocity field 
of each fluid remains potential, the total velocity field is 
no more potential. Indeed, we expect that it develops 
a rotational part due to the presence of the dissipative 
term a ij l22l. 



B. Adiabatic and isodensity modes 

As we did for the single-fluid case, let us study the lin- 
ear evolution of the multi-fluid system by dropping the 
right-hand side of eqs. (48 ) and (49 ). Since at linear order 
there is no anisotropic stress <7 U , which is second order 



in the velocities, the two linear solutions ( f ) and ( 13 ) 



found in the single-fluid case are expected to be also solu- 
tions of the linear multi-fluid system. This corresponds 
to the case where all the fluids start comoving and, as 
they all follow geodesic motion, remain comoving during 
their entire evolution. Analogously to the jargon adopted 
in the physics of the early universe, these solutions cor- 
respond to the so-called growing and decaying adiabatic 
modes. Note that if only these two modes are initially 
excited, the right-hand side of eq. ( 52 1 vanishes and the 



total fluid is indistinguishable from a pure dark matter 
fluid. 

However, the presence of multiple components gives 
birth also to isocurvature or rather, given our scales of in- 
terest, isodensity modes. To examine their properties, let 
us turn to the equations describing the multi-component 
system, eqs. (48) and (49). In this case it is convenient 



to introduce a multiplet ^f a which generalizes the duplet 
defined in eq. (12), i.e. [T4"] 



(j 1) e 1 ,j 2 ,e 2 ,...) J 



(53) 



where 9 Q = —8 a /f+(t). Thus, for N components $> a 
has 2N elements. Equations (48) and (49 1 can then be 



rewritten as eq. (15) where in this case the matrix ele- 



ments of fl a b are given by 

^(2p-l) (2p) 



n 



(2p) (2p) 



= -1 , 

2 7* 



- 1 



n 



(2p) (2,-1) 



+ 

3 flm n 



(54) 



for any integers p and q running from 1 to N and where all 
the other elements of fl a b vanish. The non-zero elements 
of the coupling matrix j a b c are 



7(2p-i) (2p-i) (2p)( k , ki,k 2 ) = 
7(2p-i)(2p)(2p-i)(k, ki,k 2 ) = 



a(k 2 ,ki) 
2 

a(k 1 ,k 2 ) 



(55) 



7(2 P )(2p) (2p)(k,ki,k 2 ) = /3(ki,k 2 ) 



for any integer p. Note that there are no explicit cou- 
plings between different species in the 7 b c -matrices. 

The isodensity modes are obtained under the con- 
straint that the total density contrast vanishes, i.e. 5 = 0. 
Since the evolution equations decouple under this con- 
straint, the time dependence of these modes can be easily 
inferred. One solution is given by 



©i is) (ry) ocexp 



drf 



3 O rn 

27T 



i 



with 



/a©i iS) = 



(56) 



(57) 



which automatically ensures that ^ Q f a Sa — 0. Note 
that as O m //f departs little from the value taken in an 
EdS cosmology, i.e. Q m //? = 1, the isodensity modes 
are expected to depart very weakly from 

(t?) a exp(-77/2) , e } (V) = -2e< is) (ry) . (58) 

A second set of isodensity modes is given by 

9(f) (rj) = , Si ci) (ry) = Constant , (59) 

again with 

E/«^ ci) = - ( 6Q ) 

a 

To be specific, let us concentrate now on the case of 
two fluids and assume an EdS background. In this case 
the growing and decaying solutions are proportional, re- 
spectively, to 



x«= (1,1,1, if , 
4-)=(l,-3/2,l,-3/2f 



(61) 
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Moreover, the isodensity modes are proportional to 



,(is) 



(-2/ 2 ,/ 2 ,2/i,-/i) J 



i d) = (/j.o.-/i.of 



(62) 



We are then in position to write down the linear prop- 
agator gabiViVo) satisfying eqs. (21) with (|22|). For two 



fluids and an EdS background it reads [7] 



/3/i 2/i 3/ 2 2/ 2 

3/x 2/x 3/ 2 2/ 2 

3A 2/ x 3/ 2 2/ 2 

\3/i 2A 3/ 2 2/ 2 

/0 -2/ 2 

h o 

2/i 

\0 -/a 



e 2 



/ 2/i ~2fr 2/ 2 -2/ 2 

-3/i 3/i -3/ 2 3/ 2 

2A -2/ x 2/ 2 -2/ 2 

V-3A 3 A -3/ 2 3/ 2 



/ h 2/ 2 -f 2 -2/ 2 



-/i -2/i A 2/i 

\ 



(63) 



In the following we explore how this propagator is 
changed by the coupling with the long-wavelength modes 
in the eikonal approximation. 



C. Resummation of the propagator with the 
eikonal approximation 

Let us study the resummed propagator in the presence 
of more than one fluid. For simplicity, we will restrict 
the study to the two-fluid case and an EdS background. 

The eikonal equation, eq. ( 28 ) with ( 29 1 , also holds in 
the multi-fluid case. 



Note that, because of cq. (59), the constant isodensity 



However, in this case is given 
by a sum of adiabatic contributions, for which the fluid 
displacements are the same, and isodensity contributions, 
for which their weighted sum vanishes, i.e., 

E ab (k, v ) = zW(k, V )6 ab + 3<?(M) , (64) 

where takes the form 



-(is) _ „( is ) , 
^ab ~ " 



(h o 
o h 

o o -A o 
V o o o -a 



(65) 



If we assume 2 a b to be in the linear regime, then 

~( ad )(k,,?)EE f d 3 q^(e(+)(q,r ? ) + e(-)(q,r / )) , 

(66) 

where 0( + ) and 1 - ' are, respectively, the growing and 
decaying adiabatic modes of the long-wavelength dis- 
placement field. The isodensity contribution 5 ^fj c on- 
tains the decaying isodensity mode given in eq. ( 58 1 , so 
that it reads 



n(is) = 



1 



^3 k q 



ef 8) = 



i 

A js 



d 3 q ^e^ 



mode does not contribute to 2 a &. 

We are now interested in computing the resummed 
propagator in the eikonal approximation under the mod- 
ulation of the long- wavelength modes in eq. ( 64 ) . As the 
adiabatic modes in S a 6 are proportional to the identity, 
their effect can be incorporated in exactly the same man- 
ner as in the single-fluid case. The adiabatic modes will 
contribute to the resummed propagator by a multiplica- 
tive factor of the exponential of the adiabatic displace- 
ment field, as in eq. ( 33 ) , 



x exp ^ dr/E^k, 77')) 



(68) 



Note again that, as in the single-fluid case, the soft adia- 
batic modes induce a phase change but do not affect the 
amplitude of the small-scale modes. 

Including the isodensity mode in the resummed prop- 
agator proved difficult. We have not been able to find 
a closed analytic form for it. Thus, we have to rely ei- 
ther on numerical studies or on perturbative calculations. 
As an example, in Fig. [3] we show the effect of the soft 
isodensity mode on the small-scale modes, by plotting 
the evolution of the resummed CDM and baryon density 
modes with s( ad ) = 0, i.e., 



S c (r), S< is >) = £ la (fc; 77, »&,; S< ad > = 0)* a (r? in ) 
6 h ( V ,E^) = £ 3 a(fc; 77, r?i„;S( ad ) = 0)* o (rfa0 



(69) 



(67) 



normalized to the growing mode e v . Initial conditions 
are chosen such that ^>a(Vin) = (1, 1, 0.15, 0.15) T and we 
have taken |S( is )(fc,?y in )| = 25. At early time the CDM 
and baryon density mode grow slower than e 71 (upper 
panel) and since |S( 1S )| 1 the phases of the two modes 
rapidly evolve (lower panel). At late time the phases 
are fixed and the density modes evolve according to the 
standard adiabatic growing mode. Let us study these 
two limiting behaviors. 
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,5 - 



Exp(v c ?7) 



l<S»l 
Exp(v,,r]) 



where 



in 




FIG. 3: Evolution of the amplitudes (upper panel) and phases 
(lower panel) of the resummed density modes with ~' ad ' = 0, 
i.e. <5 c (f?,H (is) ) and S b {rj,E (is) ) as defined in eq. (p) for CDM 
(thick solid line) and baryons (thick dashed line). The initial 
conditions are chosen such that ^a(?7m) = (1, 1, 0.15, 0.15) T 
and |E' ls '(?7i n )| = 25. Thin lines represent the analytic solu- 
tions at early times of eq. (75 1. 



1. Early-time behavior 

As S( 1S ) is a decaying mode, it can become arbitrarily 
large at early time. Let us consider a mode k for which 
initially |S( 1S )| 1. This means that k • v$, i.e. the dis- 
placement field of the soft modes along k, is much larger 
than the Hubble flow. In other words, the time scale of 
the motion of the large-scale modes is much shorter than 
the time-scale of growth of the small-scale ones, set by 
the Hubble time. 

We can grasp the nature of the early-time evolution by 
making the following change of variable, 



drf E aa ( V ') 



(70) 



where there is no summation over a in S aa . In this case 



the first line (i.e. the large-A: part) of eq. (28) can be 
rewritten in terms of ty a as 



d = 



* Q (ry) + 0,6*6(7,) = E ab ( V )^! b (r]) , (71) 



a, 



and 



with 



/ -1 

-3/i/2 1/2 

-1 

V -3/ 2 /2 1/2 





3/ 2 e~^/2 



3/1^*72 



<p(n) = —i df/S (is) (?/) 



(72) 



(73) 



(74) 



The two fluids are only coupled through s a 6. However, 
£( 1S ) j n e q (74 1 jg purely imaginary: At early time the 



coupling term contributes to a rapidly changing phase 
ip. When the time scale of these oscillations is much 
shorter than that of structure growth, this force term 
can effectively be neglected and the different species de- 
couple. Indeed, the velocity difference (in the direction 
along k) between the coherent flows of the two species 
is large enough that the short modes of one fluid do not 
gravitationally see those of the other fluid. 



The system we are left with is given by eq. ( 71 ) with 
vanishing right-hand side. For an EdS background the 
solution of this equation is given by [TS] 



5 a oc exp(i/^ ±) 7?) 



with 



24/ Q 



(75) 



(76) 



The growing solutions in eq. ( 75 ) explain the early-time 
evolution shown in Fig. [3] At early time, both CDM 
and baryons grow slower than e v and their phases are 
dominated by their respective large-scale isocurvature 
displacement fields, as accounted for by the change of 
variable (70 1. 



2. Late-time behavior 

As the isodensity mode decays, one expects to 
recover at late time the single-fluid propagator 

dabiVi Va) ex P ( J^ d?7'S( ad ' (yi'fj . More precisely, we can 
compute how the propagator £ a 6 deviates from the adia- 

(is) 

batic one with a perturbative analysis. Indeed, since S^ b ; 
becomes small at late time, one can compute £ a 6(k; rj, r/o) 
perturbatively in S( IS ). 

Solving the first line of eq. (28) at first order in S( IS ) 
yields, 
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£ ad (k, v , m )^cxp[ / dr/S( ad )(k,r/) 

»7o 



9ad(v,Vo)+ I dr ) 'g ab (r),T] / )E^ ) Ck,r] / )g cd (T] / 1 T] ) 



'Jo 



(77) 



By plugging in this equation the expressions of the linear propagator g ab from eq. (63) and of ^ ab from eq. (65), and 
integrating in time yields 



exp 



d7?'S( ad )(k,7 ? ') 



'Jo 



(78) 



with 



C ab {v,Vo) = fifa 



2e r '-no 



h(l-rio) 




( 12 8 12 77 8 £ \ 

3 2 34 2^ 

-124 -84 -12 -8 

/2 J2 

V -34 -24 -3 -2 / 

v . J2 J2 



+ 2 



/ f-3 2^-3 1-34 2-34\ 


- 1 34 - 2 3 

J 2 11 



h. — i n ^ _ oh 



3-24 

J2 



+ e ~h(v-vo) 



+ 2e 



-('J-'Jo) 



/0 2f-3 3- 2f \ 
3-4 f-3 
24 -3 3-2^ 



\0 3 



g-K^-^o) 



2 e - 2 (')- r 'o) 



/2 . . f- 

/2 U /2 
/ ^ 1 ^ ft \ 

i -l 4 

h h 
LX — IX 1 _1 

/2 /2 
Jl Jl 1 



3 / 



1 / 



/4-24 8-84 44-2 84 -8\ 

ji ji ji ji 

-2 2^-4 1-4 4-2^ 

2-44 8-84 24-4 84-S 

H 11 J2 J2 

4 - 1 24 -4 2 4-24 / 

J2 / 



/2 



/2 




(79) 



Note that this result is written in terms of S^ ls ' taken 
at the initial time iji n = 0, so that the time dependence 
of S( 1S ) is included in eq. (78) and in the square bracket 



of eq. (79). As expected by eq. (65), the corrections to 



the propagator due to the isodensity mode are invariant 
under exchange of fa O —fa and 1,2 <H> 3, 4 in the matrix 
indices a,b. Furthermore, CabiViV) — 0- 

The final expression of the nonlinear propa- 
gator is obtained after the ensemble average of 

S< is )(0) exp ( dr/S( ad )(r/)) has been taken. We recall 

here that the different modes that enter in S a j, are not 
statistically independent. The ensemble average can be 
written as (see appendix [A|, 



^)(0)e^ d ^ d) ^)^E^exp(E 



\q=2 



where c q is the g-order cumulant of the adiabatic modes 



and X\,p—i is a p-order cross-cumulant defined as 

p-i\ 



(81) 



The explicit values of such coefficients depend on the 
precise model. For Gaussian initial conditions only c-i 



and Xi 7 i are non-zero. Then, eq. (80) can be rewritten as 



S (is) (0) exp ' d?/S (ad) (r/)^ = Xl>1 exp 



where 



C2 



dj/dr/^S^r/)-^ (?/')) , 



(JO 



(82) 



(83) 



= / d7/(S( is >(0)~( ad )(7/)) . 



'Jo 



At late time ^,( ad ) is dominated by the growing mode. 
Thus, we can express it as on the right-hand side of 
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eq. (37) and we can use eq. (39) for For x± } i we 
find 



Finally, the ensemble average in eq. (82) can be written 

as 



= -fcV x (e"-e" ) 



(84) 



where <r 2 is the cross-correlation between the initial iso- 
dcnsity and the adiabatic modes, 



2 _ 1 [ J 3„ C 'm(g) 

ct x = o / d q— ^~ 



(85) 



with Cj n defined by 

y is >(q,0)e( ad )(q',0)\ = 5 D (q + q')C in ( (? ) . (86) 



S (i S ) e /," dVH(-»(r ) ')\__ t2 ^ 2 uv _ 



(87) 

so that the nonlinear propagator reads, at first order, 



G ab (k; 77,770) « [ 5a6 ('7,%)-fc 2 ( 7 2 x ( e ''-e f ' n )e-'' / 2 C ab (7 7 ,7; ) 

I 



o -feV2(e"-e"0) 2 /2 



It is possible to compute the nonlinear propagator at 
higher orders in S( IS -*. In particular, in appendix |b| wc 
derive a recurrence formula for the most growing mode 
of the resummed propagator, to any order in "E.( m \ We 
are now in the position to illustrate the effect discussed in 
this section in a practical case, i.e. the mixture of baryons 
and cold dark matter after decoupling. 



IV. CDM AND BARYONS AFTER 
DECOUPLING 

As an application, in this section we consider the case 
of baryons and CDM particles just after decoupling. This 
situation is illustrative of the concepts that we introduced 
in this paper. Here we focus on the behavior of the propa- 
gators on scales which are interesting for PT calculations, 
i.e. k < 1/iMpc -1 . We will see that for such statistical 
objects and such scales the impact of isodensity modes 
is very small. We leave the calculations of power spectra 
for further studies. 

The first step of our analysis is to properly identify the 
isodensity modes after recombination. We will assume 
that the primordial (i.e. before horizon crossing) large- 
scale perturbations are strictly adiabatic. In this case 
each fluid component is proportional to the same random 
field, for instance the primordial curvature perturbation 
C(k). We can assume that at the initial time 77j n = the 
different fluid variables are in the linear regime. Then, 
they can be written in terms of the initial linear transfer 
functions T a (k,0) as 



*„(k,0) =T o (fc,0)C(k) ■ 



(89) 



We will use CAMB [16] to generate the CDM and baryon 
initial transfer functions, assuming the following cosmo- 
logical parameters: fl c — 0.233, 51b = 0.0461, ft\ — 



0.721, h = 0.700, n s = 0.96, A c = 2.46 • 10~ 9 and mass- 
less neutrino species. 

A remark is in order here. Cosmological fluctuations, 
such as those described by the CAMB code, obey lin- 
ear general relativistic equations. On large scales, i.e. on 
scales comparable with the Hubble radius, these equa- 
tions may considerably deviate from the Newtonian equa- 
tions used in RPT, also at the linear level. Thus, one may 
worry that the transfer functions generated by CAMB 
will be affected by these deviations, which are gauge de- 
pendent. However, as shown in appendix [Cj for a set of 
pressureless fluids there exists a choice of variables for 
which at linear order the relativistic equations exactly 
reduce to the Newtonian equations. For the density con- 
trasts of cold dark matter and baryons, this choice cor- 
responds to take the energy density perturbations in a 
gauge comoving to the total fluid. For the velocity diver- 
gences this corresponds to take them in the longitudinal 
gauge. In the limit where we can neglect radiation en- 
ergy and momentum, the dynamics of these variables is 
well described by the Newtonian equations even on super- 
Hubble scales. 

Finally, note that even though J7a ^ 0, we will use the 
linear propagator derived in sec. |III| in a EdS universe. 
Indeed, as explained in [3] most of the cosmological de- 
pendence is encoded in the linear growth function D + 
and using the propagators derived for an EdS universe is 
a very good approximation. 



A. The linear modes after decoupling 

In Fig. g] we show the transfer functions for the dif- 
ferent fluid variables normalized to the transfer function 
of the total matter perturbation S m . We choose redshift 
z = 900 as initial time 77j n = 0. At this redshift the energy 
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0.001 0.005 0.010 0.050 0.100 0.500 1.000 

k[h Mpc -1 ] 

FIG. 4: Shape and amplitude of the transfer functions at 
z = 900. The transfer functions are plotted in units of the to- 
tal density transfer function. From top to bottom we have the 
CDM density transfer function (continuous line), the CDM 
velocity transfer function (dotted line), the baryon density 
transfer function (dashed line) and the baryon velocity trans- 
fer function (dotted-dashed line). On super-horizon scales 
they are all approximately equal, denoting that /+ ~ 1. One 
can observe that at this high redshift the baryon transfer func- 
tions are highly suppressed. 



density and momentum density of the radiation are still 
important (of the order of 20% percents). However, we 
will neglect their contributions in our treatment. More- 
over, since on super-horizon scales the transfer functions 
are all approximately equal, then / + ~ 1 and we will take 
Q = — 9 a . Note that, contrarily to what done in [7], one 
cannot consistently assume that the density and velocity 
transfer functions are the same. 

The linear evolution of each mode can be constructed 
by applying the linear propagator g a b given in eq. ( 63 1 
In particular 



Hub 



,(-) 



(is), 



9 a b (V,Vo), 9ab (v,Vo) and 
9^ab\ r ]i r ]o) are the growing and decaying adiabatic, and 
the decaying and constant isodensity time-dependent 
projectors, given respectively by the first, second, third 
and fourth term on the right-hand side of eq. (63). In 



terms of these projectors one can define the transfer func- 
tion for each mode as 



(90) 



T(+>(M) = (/c,0,/ b ,0) g^\v,0) T b (k,0) , 
T^ik^) = (f c ,0J h ,Q) g^i^O) T b (k,0) , 
T<M(k, rf) = (0, 1, 0, -1) g<£> (77, 0) T b (k, 0) , 

T^{k,n) = (1,0,-1,0) g^(v,o)n(k,o) . 

These definitions have been chosen in such a way that 

T a = T^ui +) + T<->u(-> + T^u^ + T^ui ci) . (91) 

These quantities are shown on Fig. [5] where we plot 
the amplitude of the transfer functions T^~\ T^ ls ' and 
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0.050 0.100 
k[h Mpc -1 ] 



0.500 1.000 



FIG. 5: The transfer functions at z — 900, normalized to the 
adiabatic growing mode. From bottom to top, the adiabatic 
decaying mode (dashed line), the isodensity decaying mode 
(continuous line) and the isodensity constant mode (dotted 
line). 



T( C1 ) at initial time, normalized to the amplitude of T^ + \ 
Note that from these results one can compute the r.m.s. 
of H( 1S ) that appeared in the previous section. One finds 
that 



;(is)' 



1/2 



= 8.6-10" 



h Mpc" 



(92) 



at redshift z = 900 

-l 



showing that for our scales of inter- 
est, k < 1 /iMpc" 1 , the effects of the isocurvature modes 
can only be small. The explicit dependence of the propa- 
gators on the isodensity modes is shown in the following. 



B. The nonlinear propagators 

In the presence of the decaying isodensity mode, the 
resummed propagator is not any more proportional to 
the free field propagator. The effect of the isodensity 
mode on the resummed propagator is modulated by the 
matrix C a b(r],r]o) in eq. (79). To show this modulation, 



let us define the quantities 

RsMv) = C la (ri,0)T a (k,0)/T 5c (k,ri) , 
R 5h (k,v) = C 3a (r],0)T a (k,0)/T 5b (k,r)) , 



where 



Ts c {k, rf) 
T s b (k,r]) 



9ia(v,0)T a (k,0) , 
g 3 a(v,0)Ta(k,0) . 



(93) 



(94) 



In Fig. [6] we have plotted these quantities as a function 
of scale and for different redshifts z — 500, 200, 40, 10, 3, 
corresponding to D + = 1.99,5.26,25.71,95.47,258.43. In 
Fig. [7] we have plotted R§ c and Rs b as a function of red- 
shift and for different scales k — 0.001, 0.01, 0.1, 1 in units 
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FIG. 6: The quantities R Sc = Ci a {z,z in = 900)T a (k,z in = 
900)/Ts c (k, z) (upper pannel) and Rg b = C3 a (z,Zi n = 
900)T a (fc, Zi n = 900) /Ts b (k, z) (lower pannel) as a function 
of scale at different redshifts. 



of h Mpc -1 . At small redshift (large 77) Rg c and R$ b are 
dominated by the most growing mode of the matrix C a f,, 
i.e. the first term in eq. (79) which grows as e v , so that 



they are independent of redshift. At higher redshift the 
decaying modes in the matrix C a b become important and 
for z — 900, corresponding to the initial time r)- ln = 0, Rg c 
and R$ b go to zero. Note that R$ b becomes infinite twice 
around z ~ 700. This is because at early times, right af- 
ter recombination, the linear baryon density contrast is 
positive and its decaying isodensity mode dominates over 
the growing adiabatic mode. Later on it reaches a nega- 
tive minimum where the growing adiabatic mode starts 
dominating. Thus Tg b (Jt, rf) crosses zero twice. 

The entire effect of the isodensity mode on the prop- 
agator is represented by the second term in the square 
bracket in eq. 



(95) 



which for 770 = ?7i n = is 



k 2 r x a 2 d (D + (r,)-l)C ab (r,,0) 



Here the parameter r x is the ratio of the isodensity- 
adiabatic displacement cross-correlation a\ to the vari- 
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FIG. 7: The quantities Rs c (upper pannel) and Rs b (lower 
pannel) as a function of redshift at different scales. 



ance of the adiabatic displacement field cr^, 

_ a\ _ Jd\nqT( is \q,0)T(+\q,0)q n °- 3 
Tx = of ~ ' Jdlnq [T(+){q,0)} 2 q n *- 3 



(96) 



(In the second equality we have neglected the adiabatic 
decaying mode.) At z = 900 this is r x ~ 0.85. For 
kadD + ^> 1 the nonlinear propagator goes quickly to 
zero, due to the exponential damping in eq. ( |88[ ). Thus, 
the key quantity responsible for suppressing the effect at 
low redshift is actually the time dependence in eq. (95) 
given by D + (z) — 1. For scales k < fed 



the corrective term ( 95 ) is always found to be extremely 
small. 

In order to be more quantitative, in Fig. [8] we show the 
effect of the isodensity mode on the nonlinear propagator 
by plotting G$ c / G^ d ^ — 1 and Gs b / — 1 as a function 
of redshift and at fixed scale fed, where 



GsMv) = G la (k;r,,0)T a (k,0) , 
G Sb (k,r)) = G 3 a(k;r),0)T a (k,0) , 



(97) 



and G 



f d) (fc,z) and Gf d \k,z) are the same quantities 



in the adiabatic case - i.e. for S( 1S ) = 0. We have used 
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FIG. 8: Effect of the isodensity mode on the nonlinear prop- 
agators normalized to the adiabatic nonlinear propagators, 



i(ad) 



SGsJG^ = GsJGfV -1 and 5G s JGf^ = G Sh /G 5 , 

where G s and G^ d ' are the CDM and baryon nonlinear 
propagators in absence of isodensity mode, computed at fixed 
scale fed, as a function of redshifts. Note that the effect on 
the CDM propagator is plotted with the sign changed. 



that the value of the variance of the displacement field 
is fi d ~ 9.2 x 10 -3 /i -1 Mpc at z = 900. For z < 50 we 
find that the effect is less that ~ 1% and for z < 9 less 
than ~ l%o. Note that the effects are of different signs 
between CDM and baryons. 

Another example where these effects could be signifi- 
cant is when the isodensity modes are set at much lower 
rcdshift. This is potentially the case for massive neutri- 
nos. However, massive neutrinos cannot be fully consid- 
ered as nonrelativistic particles during their cosmologi- 
cal history as their behavior is determined by a whole 
set of extra modes, such as pressure fluctuations and 
anisotropic stresses. We leave the study of this special 
case for the future. 



V. CONCLUSIONS 

The eikonal approximation provides an efficient for- 
malism within which exact resummation in the high-fc 
limit can be performed explicitly or numerically. We were 
able to recover the standard results obtained in [^1 El El E] 
concerning the nonlinear 2-point and multi-point propa- 
gators describing the gravitational instabilities of a single 
pressureless fluid. In particular, the propagators are cor- 
rected by an exponential cut-off whose scale is fixed by 
the amplitude of the displacement field along the wave- 
mode k. We have shown this irrespectively of the growth 
rate of the displacement field and of whether it follows 
the linear regime, thus extending the standard results 
previously quoted. 

Note that this formalism is based on a mode separa- 
tion between large-scales and small-scales. Indeed, we 



ft, 
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FIG. 9: Mode contribution per log interval to the variance of 
the displacement field from the adiabatic modes, eq. (401, at 
z = 0. 



have assumed that the large-scale modes with momen- 
tum q, collected in the random variable S, are much 
smaller that the small-scale modes k. In Fig. [9] we show 
the contribution from adiabatic modes to the variance of 
the displacement field, per logarithmic interval. As 
one can see, most of the contribution comes from modes 



with q < 0.1 h Mpc" 
q w 0.1 ~ 1 h Mpc 



but that of smaller modes, with 



1 is not negligible. This suggests 
0.3 h Mpc -1 , a better description of 



that, for k ps 0.1 
the damping could be obtained by setting a UV cutoff 
for q in eq. ( 40 ) . 



We have then extended the eikonal approximation to 
multiple pressureless fluids. In this case one can identify 
two types of modes: Two adiabatic modes and two iso- 
density modes per added species. Isodensity modes are 
responsible for new effects. Indeed, their large-scale flow 
changes the phase but also (unlike the adiabatic modes) 
the amplitude of small scales. Thus, the growth of struc- 
ture and consequently the amplitude of propagators and 
spectra are affected in a more complex way than in the 
purely adiabatic case. In this paper we focus our results 
on the propagators, leaving the study of power spectra 
for future work. 

In contrast to the single-fluid case, where the effect of 
large-scale adiabatic modes can be taken into account an- 
alytically, for the isodensity modes we have not been able 
to find an analytic form for the resummed propagator. In 
this case, one should rely on a numerical or a perturbative 
approach. The latter is sufficient when one considers the 
case of CDM-baryon mixing. For this example, we found 
that the impact of isodensity modes on the propagators 
is very small at low redshift and for scales of interest for 
standard PT, i.e. for k < lh Mpc -1 . However, there 
might be cases where the impact of large-scale modes is 
more significant, for instance when the scale of interest 
are close to the non-linear regime at the time the isoden- 
sity modes are set in. This is expected to be the case 
for massive neutrinos. Although we did not address this 
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case explicitly, we stress that the eikonal method can be 
used irrespectively of the field content of the system. We 
leave the case of massive neutrinos for further studies. 

Acknowledgements: FV wishes to thank Toni Ri- 
otto and Ravi Sheth for interesting discussions. 



Appendix A: Computation of moments 

We want to compute the value of 

ffn = <H n exp(Z>)>, (Al) 

where S and V are two random variables whose statis- 
tical properties are entirely characterized by their joint 
cumulants, 



(A2) 



It is convenient to introduce the auxiliary function 
exp(2? + AS) and to notice that 



(A3) 



The ensemble average that appears in this expression 
can be written in terms of the cumulants generating func- 
tion of V + AS as 



A=0 



which can be re-written as 



g n = cxp ( V ) ^— exp ( V ^y\ q 



A=0 



where 



oo 1 



p=q 



(p- «)!•*"-«• 



(A4) 
(A5) 

(A6) 



g n can then be formally expressed in terms of X q as 



5o = cxp V , 
9i = Xi g , 



92 



(Xf+X 2 ) g , 



53 = (^ 1 3 + 3A 1 A 2 +A 3 ) 50 



(A7) 

(A8) 
(A9) 
(A10) 



Note that for a unit 50 the relation between X q and g n 
is exactly the one relating cumulants of order q with mo- 
ments of order n. Thus, in general this relation is ob- 
tained by the Arbogast-Faa di Bruno formulae. 

These relations greatly simplify in the case of Gaussian 
initial conditions, assuming - without loss of generality 
- that xo,o = x i.o = x o,i = 0. Indeed, in this case only 
X\ and X 2 are non-zero, 



Xi = xn , X 2 = x 2 fi , 



so that 



5o = cxp 



V 2! J ' 



91 = x i.i go , 

92 = +x 2fi ) g , 

93 = ( x i.i + 3:ei,iX2,o) 50 



(All) 



(A12) 

(A13) 
(A14) 
(A15) 



Appendix B: Higher-order time dependence 



It is possible to compute the nonlinear propagator at 
higher orders in S( 1S ) using the expansion 



where 



£ af (k, V , m ) = g af ( V , m )exp^ dv'Z^ik,^ +C { a ) ) (k, V , m )+^ a 2 f ) (k, V , m ) + ... , (Bl) 



ee exp f [" dr/S^k,,/)) r^'g^M^^gcfiri'M , 

\Jrto / Jno 



no 

V rn' 



no J no 



eg^cxp / dr/S^k,,/) / drf d^^^OS^^k^Offed^^'O^^k,^)^/^',^), 



(B2) 



In particular, here we derive an explicit expression for formula, 
the most growing solution to the, in terms of a recurrence 
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The key point is to show that the fastest growing mode (77,770) oc e A '' 7 ~' ?0 ^ and 3*^ (?y) oc e~^ v {r)i n — 0), its 

of the amplitude of the corrections to the resummed prop- time dependence is proportional to n nested integrals, 
agator goes as the adiabatic growing mode, i.e. oc e^"* 70 , given by 
to all orders in S( 1S ). Let us consider the amplitude of 
the n th -order correction. From eq. (B2) and using that 



j(n) 



-,<n 



>V2 



dr)„ 



l(»Jn- 



(B3) 



where each A^ can take the values {1, 0, —1/2, —3/2}. As 
we are only interested in the fastest growing mode, we 
take Ao = 1. 
We define 



= Xj - Xi - 



3 - ? 



(B4) 



satisfying aij + ctjk — ot-ik- It is important to note that 
a oi < for i =/= 0. Then, we can rewrite the integrals 
above as 

rv rm fhi-i 

C I dyyi e aoim / dn 2 e ai2ri 2 / dr?„ e "- 1 -""" 

Jt) Jr/o Jr/o 

. (B5) 

with C = e XaV e XnTI ° . By performing an integration by 
parts on rji, we find the boundary term 



C 



/ dr, 2 e" 12 " 2 • • • / dr?„ e"- 1 """ 

«oi J m> J V0 



and the remaining integrals 

rv potoini 

-C dr)i e ai2m ••■ 

J Vo "01 



(B6) 



Vn-l 



dr) n e a -^" . (B7) 



'/n 



The lower bound for the boundary term obviously van- 
ishes. Therefore, the boundary term contains an e a ° lV 
factor making it subleading with respect to the remain- 
ing integrals, which for large 77 become constant. To see 
this, we perform (n — 2) more integration by parts, ev- 
ery time dropping the boundary term for the very same 
reason, until we are left with 



(-1)* 



C 



dn n e Q °""" 



The leading term of this integral is independent of 77, 



(B8) 



(-irn 



c 



,ao„J7o 



(B9) 



and by replacing C by its definition, we have 

n 

j(n) = TT e -%Vo e V-r)o_ (BIO) 



«i0 



Now that we have studied the time dependence, we 
can reintroduce the time independent matrices h a bgl^ in 



( B2 ) and sum over A^ in each of the integrals of eq. ( B3 1 



Let us define the time independent matrices 



AO) = (i) 

^ab — ilab ' 



.4 



Kc9 c 



W _ V 

ab ^ l + i/2-X 

A=l,0,-l/2,-3/2 ' 



% > 1 



(Bll) 



where h ab is defined in eq. ([65]), gj£ = gil\rj , ry ) are 



the time independent projectors on the right-hand side 
of eq. (63l, and the sum runs over A = 1,0, —1/2, —3/2 
corresponding to the 4 linear modes +, ci, is and — . The 



denominator (1 + i/2 — A) is exactly the of eq. ( B10[ ). 
It is easy to verify that the most growing solution for 
in eq. (Bl I is then given by 



C\^V,Vo) = [S( is )(k,0)e-W 2 j [A(°>n A( 

i=l 

x e"-"° exp ( f dr7'S( ad )(k, rf) 

Viio 



ab 



(B12) 



This equation formally generalizes the expression for the 
most growing mode contained in eqs. (78 1 and (791 to 



any order. Note that the time dependence is the same as 
the linear growing mode, i.e. oc e n . Ensemble averages of 
these quantities can be taken using the equations given 
in appendix |A") 



Appendix C: Treating super-horizon scales 

Let us consider a linearly perturbed FLRW metric in 
longitudinal gauge with only scalar perturbations. In the 
absence of anisotropic stress the two metric potentials are 
identical and the metric simply reads 

ds 2 = a 2 (r) [-(l + 20 lon )dr 2 + (l-20i on )dx 2 ] , (CI) 

where r is the conformal time defined by dr = dt/a(t). 
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Combining the 00 and the Oi components of Einstein's 
equations gives, in Fourier space, [17] 



fc 2 ^lon = -^% 2 (fton 



3H 2 e iOD /k 2 ) 



(C2) 



where S\ on and 9\ on are the total density contrast and di- 
mensionless velocity divergence [23], respectively, in lon- 
gitudinal gauge and H = dlna/dr is the conformal Hub- 
ble rate. Moreover, in this gauge, the continuity and 
Euler equations for a single-fluid read, at linear order, 



si 



Ion 



-Wi on + 3< n , (C3) 
-(H'/H + H)9 lon + k 2 ^ on /U , (C4) 



where a prime denotes the derivative with respect to con- 
formal time. 

One can check that, at the linear level, the Newtonian 
Euler equation ^ is the same as its relativistic version, 
eq. (C4), while the continuity equations ([!]) and (C3) 
differ by the term 3</>f on . Indeed, for instance, using that 
a oc t 2 in matter dominance, one can check that the 
solutions to the above equations are 



and 



<kon(k,T) =^+(fc) + (fcT)- 5 0_(fc) 



2+-(fcr) 2 j 0+ 

^(fcr)- 3 -3(fcr)- 5 
6 



(C5) 



(C6) 
(C7) 



Thus, eq. ( C7 1 correctly describes the growing and decay- 



(101, 



ing solutions of 9 in the Newtonian limit, eq. (11), with 
9 + oc a and Q- oc a -3 / 2 , even on super-Hubble scales, 
while for (5i on we recover the Newtonian case, eq. 
only in the limit fcr > 1. 

One can define a quantity which obeys the Newtonian 
continuity equation even on super-Hubbies scales. The 
comoving energy density perturbation, defined as 



^lon 

1 

_ 6 



3H 2 9 lo 



Ik 2 

((fcr) 2 0+ + (kr)- 3 ^) , 



(C8) 



does this job. Indeed, replacing <p[ on using the Oi com- 
ponents of Einstein's equation [T7j, and using eq. (C4), 
eq. (C3) reads 



5' 



-ne lQ 



(C9) 



thus reproducing the linear part of the continuity equa- 
tions in the Newtonian limit, eq. ([!). Moreover, in terms 
of this variable eq. |C2| becomes a Poisson-likc equation, 



n 2 6 c 



(CIO) 



k 4>\on = 

reproducing eq. ([3|. 

We conclude that at linear level the Newtonian equa- 
tions ([!]), |2]) and ^ describe the relativistic dynamics 
once we interpret the Newtonian potential <fi as the met- 
ric potentials in longitudinal gauge 4>\ on , the Newtonian 
density contrast 8 as the comoving energy density per- 
turbation (5 com and 9 as the dimensionless velocity diver- 
gence in longitudinal gauge, 9\ OTl . 

The case of many fluids is not very different. In this 
case one can show that for each species a the relativis- 
tic perturbation variable which satisfies the Newtonian 
continuity equation is the energy density perturbation 
comoving to the total fluid, 



S a ,lon + 3H 2 9 lon /k 2 



(Cll) 



where 9i oa is the total dimensionless velocity divergence 
in longitudinal gauge, 



Aon 



f a 9 a , 



(C12) 



Using this variable, eq. (C2) becomes a Poisson equation 



while the velocity divergence in longitudinal gauge for 
each species a, 6 aj \ on , naturally satisfies the Newtonian 
Euler equation. 

Let us connect these variables with those of CAMB in 
the case of a mixture of cold dark matter and baryons. 
CAMB uses synchronous gauge comoving with cold dark 
matter [IB]. Using the gauge transformation between 
synchronous and longitudinal gauge [17] we have that 



^a.CAMB — <5a,lon + 3V.CT , 
#a,CAMB = #a,lon ~ ffc 2 /% 



(C13) 



where a is the shear. This is related to the usual syn- 
chronous metric perturbation variables h syn and ?7syn, 
respectively the trace and traceless part of the spatial 
metric, by a — {h' syn + 6' syn )/2k 2 . Using these equa- 
tions, since in CAMB the dark matter velocity vanishes, 
#c,camb = 0, we have 



#c,lon = ok /H , 
9b,\on — #b,CAMB + vk 2 /% , 
<5ct, com = <5a,CAMB + 3V. 2 /b^b.CAMB/^ 



(C14) 
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